"Vortex-melting front" in thin superconductors with pinning 
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Magneto-optical observations of a second flux front, which occurs at the second peak in the 
magnetization of P^Si^CaC^Orr single crystals related to the known first order "vortex-lattice 
melting", are reconsidered. We show that, in thin samples, electrodynamics necessarily leads to an 
extended region in which the magnetic induction adopts a nearly constant value close to that at 
which the phase transition occurs at thermal equilibrium. In this region a dynamical phase mixture 
of vortex "solid" and "liquid" should exist. Interestingly, the observed second flux front does not 
mark the "melting" front, as it was naively interpreted earlier, but indicates the disappearance of 
the last "solid droplets". 
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I. INTRODUCTION 

One of the most fascinating features of the vortex 
structure in layered superconductors is the so-called 
"melting" of the vortex ensemble. It is clearly estab- 
lished by now that this transition is a first-order phase 
transition above which the long-range correlation along 
the vortex lines is lost .El The most thoroughly investi- 
gated material, in which the "vortex melting" occurs well 
below the second critical field, is the high-T c supercon- 
ductor Bi 2 Sr2CaCu20 a: (BSCCO). There are attempts to 
add some other phase transition lines to the B-T phase 
diagram of this material (B is the local induction and 
T the temperature). For example, the field of abrupt 
increase of the screening current with increasing B, the 
so-called second peak of the magnetization, was ascribed 
to a second order phase transition connected with the 
"melting" line at a tricritical point. However inspection 
of the available data suggests that the line of the sec- 
ond peak in the phase diagram is the continuation of the 
"melting" line into the low-temperature area where bulk 
currents do not completely relax during the experiment .El 

It is natural to assume that, in spite of the macro- 
scopic irreversibility, the vortex phases are in equilib- 
rium locally. Then the unified "melting line" including 
the second peak-line, which shifts as one curve under 
the introduction of additional disorder and/or change_n£ 
anisotropy by oxygen doping or electron irradiationjfm 
is the only first-order phase transition line. Recently 
the absence of the tricritical point at the border of 



observable pinning was elegantly demonstrated by un- 
masking the equilibrium magnetization step at the vor- 
tex melting point from tke hysteretic magnetization us- 
ing ac demagnetizationfl This picture implies differ- 
ent pinning properties of these two vortex phases with 
correspondingly different current -voltage characteristics, 
which were extracted from local Hall probe measure- 
ments of the magnetic relaxation in the vicinity of the 
second magnetization peak at different temperatures and 
starting fields.B The resulting jump of the screening cur- 
rent at the first-order phase transition is then observed 
as a second front of the penetrating magnetic flux. 

Magneto-optical observations of this second magne- 
tization front initially did not bring much new ing 
formation and did not attract enough attentionfj'0 
Such experiments were initially interpreted straightfor- 
wardly as the direct observation of the phase boundary 
which was already well established by local Hall probe 
measurementsflla However, after more careful considera- 
tion it became clear that the suggested "naive" interpre- 
tation is not compatible with the (macroscopic) electro- 
dynamics of thin samples. Namely, in thin superconduc- 
tors a sharp boundary between two different screening 
currents should produce the same logarithmic peak in the 
induction as at the sample edges. However, this moving 
peak is not observed in experiments. The correct theory 
indeed yields that a transition region with finite width 
and smooth current distribution should be expected, as 
it is known from the exact solutions for the penetration 
of perpendicular flux into thin samples, see e.g. Rcf. |. 
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We propose here a model which predicts that this tran- 
sition area in thin samples should inevitably contain a 
mixture of both phases. 

A comprehensive description of the peculiarities of 
magnetic flux dynamics became particularly important 
when new observations of metastable vortex phases un- 
der fast remagnetization were reported.E30 In order to 
separate new dynamic effects from the (already compli- 
cated) geometry-dependent "regular" flux dynamics oc- 
curring when the vortex phases are in local equilibrium, 
we make a step backward in our investigations and de- 
scribe the equilibrium case more carefully. The nature 
of flux flow at the phase transition is simulated qualita- 
tively using an algorithm developed for the flux dynamics 
in thin conductors with nonlinear current-voltage law.113 



II. MAGNETO-OPTICAL OBSERVATIONS 

The penetration of a magnetic field into BSCCO single 
crystals was visualized by means of the magneto-optical 
imaging technique using a ferrimagjaetic garnet indica- 
tor film with in-plane anisotropy.£3iij In general, higher 
light intensity in the image means higher local magnetic 
induction normal to the sample surface. The images were 
recorded by a CCD camera during the continuous ramp 
of the perpendicular magnetic field with dB a /dt = 200 
G s _1 and digitalized every 250 ms by a video- frame grab- 
ber card. The camera was working at the standard TV 
frame frequency of 25 Hz, so each image of the obtained 
movie is an average over about 40 ms. The increase 
of B a during the acquisition of a single image is com- 
parable to the measurement noise. The relatively fast 
ramp of the field was used in order to observe the critical 
slope of the initial penetration, which otherwise relaxes. 
The intensity profiles measured across the digital images 
were calibrated in terms of field units using an additional 
recording under equivalent conditions but above T c . A 
separate calibration curve was obtained from the addi- 
tional "empty" images for each point on the flux profiles. 
In this way we reduce the effect of inhomogeneities in the 
illumination and in the sensitivity of the indicator across 
the image. For demonstration, an underdoped BSCCO 
crystal grown at the FOM-ALMOS center (the Nether- 
lands) by the traveling-solvent floating zone technique 
under 130 mbar O2 partial pressureEEl was specially se- 
lected for the uniformity of flux penetration and for the 
low field B sp w 180 G at which the second magnetization 
peak occurs. The crystal was cut into a rectangle with 
dimensions 1450 x 700 x 43 /xm 3 . 

The initial penetration into the zero-field-cooled sam- 
ple resembles the numerous magneto-optical observations 
of the Bean- like critical states in type-II superconductors. 
The crystal edges are well observed as a sharp change 
from the very bright field outside the sample and the grey 
interior. The flux penetrating from the crystal edges is 
visible as a grey belt (higher brightness represents higher 
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FIG. 1. Magneto-optical imaging of the flux penetration 
into a zero-field-cooled BSCCO single crystal at T = 28.8 K 
in continuously ramped B a with ramp-rate of 200 Gs _1 . 
In every image, light areas correspond to high induction, 
while dark areas have low induction; however, the inten- 
sity-to-induction mapping (look-up table) in different images 
is selected separately in order to reveal specific details. Se- 
lected snapshots from the movie: (a) B a = 200 G, ordinary 
flux front below B ap . The dark area in the center is the Meiss- 
ner area with B = 0, the bright area around the crystal corre- 
sponds to the applied field B a and the grey belt at the edges 
represents the classical Bean-like penetration with a critical 
slope of B. (b) B a = 300 G: the second flux front appears 
as a new, sharp increase of the image brightness at the top 
and right sample edges, while the initial field profile persists 
in the sample center, (c) B a = 450 G: the second flux front 
above B sp . Note that the dark central part now corresponds 
to B ~ B sp , unlike the similar frame (a), where this value is 
represented by a light grey tone. The calibrated flux profiles 
for the entire image sequence, shown in Fig. 2, are taken along 
a vertical line crossing the sample at its center downwards. 

local value of the perpendicular magnetic induction) 
which squeezes the central (dark grey) Meissner area in 
which B is zero (Fig. |l| a). The gradual decrease of 
the brightness of the penetrated perimeter corresponds 
to the well-known Bean gradient of the penetrating flux. 
When the field at the edges reaches B sp , a second sharp 
flux front starts to penetrate from the crystal edges 
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(Fig. @ b)B&B When the second flux front has pene- 
trated deep into the sample, the resulting image is be- 
coming nearly identical to the initial penetration (com- 
pare Fig. [l] a and b). However, one should remember 
that the induction in the central dark area in Fig. ^ c 
approximately corresponds to that in the light grey field 
outside the sample in Fig. [I] a, where B — in the center. 



try, where it occurs at the junction of two different field 
gradients. In perpendicular geometry, however, any flux 
front penetrating into a thin superconductor is accompli^ 
nied by strong sheet currents inside the screened areanxn 



III. MODEL CALCULATIONS AND DISCUSSION 
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FIG. 2. Magnetic flux profiles across the sample shown in 
Fig.l at T = 28.8 K. The field of the first order "vortex melt- 
ing" is indicated by two dashed lines. The position of the 
sample edges is marked by vertical dotted lines with arrows 
at the abscissa. The profiles represent a sequence of images 
taken every 0.25 s, i.e., at intervals AB a — 50 G. The images 
for the profiles number 4, 6 and 9 (counting from the bottom) 
are shown in Fig.l a, b, and c, respectively. 

The whole course of flux penetration during field ramp- 
ing from zero to 500 G is visualized in Fig. ^|by B(x) pro- 
files measured across the sample width. After reaching 
a well-determined field level the increase of the induc- 
tion B pauses and B is stabilized until a new flux front 
approaches which increases B further. The plateau of 
constant induction is associated with the phase transi- 
tion at B sp . The slight inclination of this plateau, which 
is marked by two dashed lines, is due to inhomogeneity 
of the crystal. To our knowledge, such inhomogeneity is 
found in all BSCCO crystals and is inevitable (see e.g. 
Refs. ^7 18). During the subsequent decrease of the field 



a similar feature appears at the same line B sp without 
notable hysteresis in the transition field. The second flux 
front looks very similar to what the ordinary Bean-like 
flux profile would have looked like if some flux relaxation 
had occurred during a temporal halt of the field ramp. 
Another close analogy exists with respect to flux pene- 
tration into a superconductor after field cooling. That 
analogy, in particular, shows that the initial interpreta- 
tion of the second front as a phase boundary separating 
regions with very different critical currents, is not com- 
patible with the perpendicular geometry of the experi- 
ment. Such a sharp border between two currents at a 
definite field is possible only in the longitudinal geome- 



In order to understand the experimental observations, 
let us first consider the simplest possible situation of a su- 
perconducting strip occupying \x\ < a, \z\ < d/2 (d <C a) 
in which the vortex matter undergoes a first order tran- 
sition at B z = B sp . In the vortex solid below B sp , the 
critical current is zero, and in the vortex liquid above 
B sp , it has some non-zero value J c . In other words, the 
critical current discontinuously jumps from zero to J c as 
the vortex matter melts. Furthermore, let us neglect sur- 
face and geometrical barriers. The flux penetration into 
this strip is illustrated in Fig. ||. As long as the applied 
field B a || z does not exceed B sp , the strip is perfectly 
transparent to the field (flat profiles 1 and 2 in Fig. |^). 
Once the applied field exceeds B sp , the internal induction 
is stabilized at the value B sp , and further flux penetra- 
tion is screened by a non-zero J(x). The distribution of 
this J(x) and of the internal induction are described by 
the well-known solution of Ref. [)| 
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where B c — TTfioJc/d is the characteristic field, ±6 = 
a/ cosh[(i? a — B sp )/B c ] is the position of the flux front, 
and c = \J a 2 — b 2 . For Fig. [| we chose B c = \B sp . The 
situation is illustrated by the profiles 3 and 4. In the 
region in which B > B spi the current J = J c , while in 
the region with B = B sp , the strip geometry imposes the 
smooth decrease of J to zero in the interval between the 
"second" flux front and the sample center. While it is 
evident that the former region is occupied by the high 
field vortex liquid, the latter region must also consist of 
vortex liquid, since the vortex solid (with zero critical 
current) cannot sustain the non-zero J(x). Therefore, 
the observed flux front does not correspond to the front 
of phase transformation. 

The above consideration, based on the simple Bean 
model, does not account for a non-zero current in the 
solid phase, nor for the possibility of flux creep. We in- 
clude these effects using the algorithm previously devel- 
oped for studying the flux penetration into thin supercon- 
ductors with a nonlinear current -voltage law E(J, J c ) = 
Po\J/Jc\ a J where J is the sheet current lithe current den- 
sity integrated over the strip thickness) £3 We model the 
sharp jump of the current-voltage law from E~{J) = 
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E(J, J~) in the low-field low-pinning phase (J c = J~) to 
E + {J) — E(J.J+) in the high-field high-pinning phase 
(J c = J+) by introducing a step in J C {B) at the transi- 
tion field i? s 



MB) = J, 



cO 



AJ C , B 
— — — tanh — 



B. 



sp 



AB 



(1) 



with jump width AB <C S S p and jump height AJ C = 
J+ — J~ . Other models of steps in the current -voltage 
characteristics, e.g. a jump in the creep exponent cr(B), 
or different values of the parameters, give qualitatively 
similar results. 
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FIG. 3. A simplified model for the flux penetration in the 
presence of a first order vortex matter phase transition, the 
critical current is zero for B < B ap and equal to 1 for B > B sp . 
Profiles 1 to 4 correspond to B a = B c , 2B C , 3B C , and 4B C . 
The field B sp is equal to 2B C . 

The following graphs (Fig. |i| a-c) show profiles (across 
the strip) of the magnetic field B(x), screening current 
J(x), and the electric field E(x) under a constant ramp 
rate of the external field dB a /dt = 1 for a = 20, B sp = 1, 
AJ C — 0.5 and AB — 0.04, in dimensionless units where 
J+ = 1 (critical sheet current), a — 1 (half width of 
the strip), B = H, dE/dx = dB/dt (left-hand coordi- 
nate system), and po = 1. Our computed profiles of the 
induction B(x) (Fig. |i^,) well resemble the experimental 



profiles (Fig. |2|). The initial penetration at B < B sp is the 
same as was calculated analyticallyQ and numericallyEj. 
As usual, the corresponding J(x) in the flux penetrated 
area saturates at the "spatial projection" J~ at of the low 
field current -voltage characteristic. The reason for this is 
very simple: after the complete penetration of the field, B 
increases everywhere with the same constant rate. From 
the Maxwell law dE(x)/dx = dB(x)/dt it follows that 
E~(J^a,t{x)) = xdB a /dt ex x everywhere. Thus J^ t { x ) — 
J~ (x dB a /dt), where J~(E) is the inverse function of 
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FIG. 4. Computed flux penetration into a superconduct- 
ing strip with a sharp jump of the critical current from 
Jc = 0.75 to J^" = 1.25 at the vortex-lattice phase-transition 
field B S p = 1: (a) perpendicular induction B{x), (b) screening 
sheet current J(x), (c) electrical field E(x), and (d) relative 
concentration of the two phases n(x) for B a = 0.2, 0.4, . . . 
2.6. The right dashed vertical line indicates the well-known 
initial penetration front at B a = 0.2 with a sharp decrease of 
J(x) from the critical value J s ~ t w 0.7 and E(x) dropping to 
zero. To illustrate the correlation between B(x), J(x), E(x) 
and n(x) during "vortex melting", the profiles at B a = 1.4 
are marked as bold lines. B(x) clearly exhibits a second flux 
front on top of B sp = 1, starting at the left vertical dashed 
line, where J(x) sharply decreases from J s ^ t ~ 1.2 towards 
J s ~ t w 0.7 and E(x) has a sharp kink. But the corresponding 
change of n(x) extends much farther towards the center, as 
indicated by an arrow pointing from the second front to the 
left. The inset in (c) shows the geometry. All values B, J and 
E are given in the dimensionless units J c (see the text). 
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E~(J) denned below B sp . The constant critical current 
of the classical Bean model corresponds to the nearly 
constant flat part of J^ t (x) « J s ~ t . Unlike the Bean 
model, in our more realistic model the abrupt change 
of direction of the critical current at x = 0, the so- 
called discontinuity line (d-line)j£3 has a finite thickness. 
The transition from +0.8J s ~ t to +0.8J s ~ t extends over 
Ax d w 2 £-(0.8 J-)/{dB a /dt\. For the particular E(J) 
used in our simulations (Fig. |j) Ax^ « 0.02a. 

When B increases further and approaches the value 
-B sp , the field gradient decreases, and further flux pen- 
etration occurs via the penetration of a new flux front 
moving on top of the line B rs B sp . It is important to 
note that, as in the simplified case outlined above, the 
screening current J(x) does not jump sharply at the flux 
front but changes smoothly from J s ~ t (x) to the new sat- 
uration J+ t (x) w J+ t before the second front has arrived 
(Fig. ||b). Here again E + ( J+ t (x)) = xdB/dt cx x yields 
the "spatial projection" J s "^ t (x) of the high-field current- 
voltage law. This behavior is valid for any pair of func- 
tions E~ (J) and E + (J) , which play the role of saturation 
profiles of the sheet current J(x) after the complete pen- 
etration of the field, below and above B sp , correspond- 
ingly. Preceding the second front there is a region where 
the flux density is nearly constant in space and time, 
B(x,t) w B sp (Fig. |]a); in this region the slope of the 
electric field is reduced since dE/dx — dB/dt is small. In 
order to reach its final saturation shape E(x) = x dB a /dt, 
E(x) exhibits a sharp upward bend at the second flux 
front, which is clearly seen in Fig. ^c. We should men- 
tion that a non-zero AB is included here only for the 
stability of the numerical algorithm; its value only af- 
fects the rounding of sharp features in the flux profile for 
B w _B sp , but not, as we checked, the wide smooth tran- 
sition of the current. The simplified solution outlined 
above shows that the smooth transition in J exists for 
AB = as well. 

Within the framework of our assumption that only 
two phases with fixed current- voltage laws E~(J) and 
E + {J) can exist in the superconductor, the above out- 
lined smooth transition of the resulting J(x) between the 
corresponding saturation profiles can be explained only 
by a phase mixture. One can imagine various configu- 
rations of the phase domains corresponding to different 
conditions, varying from the equality of E in the two 
phases to the homogeneity of J. In general the dynamics 
of the corresponding phase mixture should be solved, but 
this complication is outside the scope of our paper. Here, 
we present the simplest and most natural model in which 
the intermediate region consists of narrow lamellar do- 
mains of alternating phases oriented parallel to the flux 
flow, i.e. perpendicular to the current, across the width 
of the strip. This model corresponds to homogeneous J 
across the phase boundaries and a modulated E. In the 
mixed region, we then have Jh qu id = Jsoiid — J{ x ) an d 
E{x) = nE+{J{x)) + (1 - n)E-{J{x)). Here, n = n(x) 
is the relative concentration of one phase; n — corre- 
sponds to the "vortex solid" and n = 1 to the "vortex 



liquid". From the calculated J(x) and E(x) we obtain 
n(x) in the form 



n(x) — 



E(x)-E-(J{x)) 
E+(J(x))-E-(J(x)) 



(2) 



(see Fig. |d). In case of the simplest domain structure of 
the phase mixture, equal parallel lamellae along the flux 
flow, the relative width of the "liquid" domains is equal 
to n(x). An example of such domains at B a = 1.4 is 
shown in Fig. |^, which is obtained by flipping the corre- 
sponding profile n(x) (outlined by the bold line in Fig. ^c) 
several times and cutting it at the right at the position 
of the second flux front, and at the left, e.g. at a value 
n(x) = 0.01 (0.01 is an arbitrary choice, the final result 
depends on it logarithmically). 




FIG. 5. Scheme of the phase lamella preceding the second 
flux front marked here by a vertical dotted line that corre- 
sponds to the left dashed line in Fig. 3. The front moves from 
the right to the left. The "liquid phase" (right, white) inter- 
penetrates the "vortex solid" (left, light grey). The vertical 
scale is not given by our theory. 

Interestingly, "melting" [i.e. n(x) > 0) starts well be- 
fore the second flux front has arrived, and the second 
flux front corresponds not to the midpoint n = i, but 
to the point where n(x) has already well saturated to 
n = 1. At this point the "melting" transition is com- 
plete. The smooth transition of the current is charac- 
terized by the fast flux flow through the weaker "solid" 
lamella (light grey), the width of which plays the role of 
a "flux valve" , while the vortices in the "liquid" lamella 
(white) are nearly immobile until the second flux front 
arrives. At x < x sp , the average electric field equals 



E(x) w E, 



sp 



Xsp dB{x,t) 



dt 



dx « const 



B sp . The vortex flow 



because B(x,t) is fixed to a value 
through the region of the phase mixture is concentrated 
in the "solid" phase. The local electric field (vortex ve- 
locity) in the solid is very high, £ so iid = E~(J(x)) = E sp , 
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because the electric field in the "liquid" phase -Eii qu id = 
E + (J(x)) is nearly zero. Consequently, the width of the 
"solid" valve adapts itself so that the necessary average 
flux flow is conserved, i.e. 1 — n sp f» E sp / -E S oiid- 

Particularly impressive is the moment when the second 
front approaches and the high second critical current J^ t 
passes through the "weak" "solid" phase. Here all neces- 
sary flux flow squeezes into the very narrow "solid neck" 
with a huge speed corresponding to E~(J^ t ). The order 
of the introduced values illiquid and -Esoiid can be read 
from the "spatial projections" of E + (J) and E~{J) in 
Fig. (see the corresponding lower and upper asymp- 
totics Jta,t{ x ) an d Jsa.t{ x ))- For this one has to take a 
horizontal projection of, for example, the cross-point of 
the bold line J{x) (at B a = 1.4) and the vertical dashed 
line, indicating the second flux front, onto J^ at (x) and 
J~ at (x). Then the first one will give x + = illiquid ~ 0.05 
and the second one will be well outside the graph at the 
right x" — -Esoiid ~ 10 3 . The relative width of the "solid 
neck" is here 1 — n sp ~2x 10 -4 . In reality, even for such 
domains as thick as the sample the "neck" width given 
by our qualitative theory is well below the vortex spacing 
of 0.3 fxm. at B sp . This means that probably only a single 
line of "solid" vortices will move before the second front 
arrives. 

The proposed dynamical intermediate structure of the 
coexisting phases should not be confused with the or- 
dinary intermediate phase that exists at any first or- 
der magnetic phase transition in the narrow field range 
AB = B + — B~ over which the jump occurs. The relative 
size of the equilibrium domains n eq in that well known 
static case is fixed by the condition that the average local 
field equals the B(x) which is determined by the entire 
current distribution, n cq B + + (1 — n cq )B~ = B{x). 

As opposed to this, in our case the relative phase ra- 
tio n(x) is determined by the different dynamics of the 
flux flow in the two different phases and by the global 
electrodynamics of the system, see the profiles in Fig. |]d 
and Eq. (2). In the ordinary static case, the region of 
the coexisting phases may be very narrow because AB 
is much smaller than the typical variation of B(x) over 
the sample; its width equals (dB / dx)^ 1 AB at the phase 
transition. In our case the width of the transition region 
containing the phase mixture is a sizable fraction of the 
sample size, while our dynamic B(x) varies even more 
due to the presence of high bulk currents. 

The question arises whether the intermediate structure 
of dynamically coexisting solid and liquid vortex phases 
also occurs in thick samples (thicker than the penetration 
depth A). In the case of the simplified model with zero 
current for B < B sp and a critical current J c for B > B sp , 
we refer to Ref . pl| in which it is shown that the shape of 
the flux front in the interior of a thick superconductor is 
very well described by the shape of the current profile in 
an infinitely thin superconductor. Using Fig. |3|, we can 
thus immediately draw the current front inside a thick 
sample. The region between this front and the sample 
surface carries a current J c and must therefore be in the 



liquid state. From the Bean model only, there is no ob- 
jection to the inner core remaining in the solid state with 
J c = 0. In that case, one would also have a coexistence of 
solid and liquid phases, but now across the sample thick- 
ness. Such a scenario would depend on the stability of 
the curved phase boundary — any phase mixture in the 
inner core up to its complete transformation to the solid 
phase also does not contradict to electrodynamics. 




FIG. 6. Hypothetical phase front between the vortex solid 
and the vortex liquid in a thick sample in the situation where 
the solid cannot sustain any current flow, and the liquid has 
a non-zero critical current. The shape of the front just corre- 
sponds to that of the current profile 3 in the lower panel of 
Fig. | 

Recently, a real-time differential treatment of magneto- 
optical images allowed to visualize the propagating vor- 
tex "melting" front through a BSCCO crystal as a-Gop. 
tour of the related equilibrium magnetization jumpOa 
The melting point appeared to be very sensitive to tiny 
inhomogcncitics of the sample defect disorder and a very 
irregular "vortex melting" with no sign of equilibrium do- 
mains is observed even in the best crystals. Our lamella 
should be more robust because they are determined by 
the critical currents, which are much higher than the tiny 
equilibrium current drop at the "vortex melting" point. 
But unfortunately, the elaborate procedure of such obser- 
vations cannot be applied to our case of irreversible mag- 
netization. Very recently, after preparation of this paper, 
a mixture of phases with different E(J) has been ob- 
served in the peak effect region in NbSe2 using a s cann ing 
Hall probe magnetometer and local ac excitation.E3c3 We 
have to mention that the phenomenon observed there, is 
static and the suggested model of the "vortex capacitor" 
differs from our dynamic model. Furthermore, we deal 
with an equilibrium transition, contrary to the metastable 
state considered in Refs. ^2] and |2^. The new observa- 
tion technique, nevertheless, may be applicable to our 
case; its magneto-optical imaging counterpart can be re- 
alized as well. Note however that in our case the picture 
is evolving rapidly (on the time scale of 10~ 2 s) and the 
irreversible currents which would screen the small ac ex- 
citation are much higher. 



IV. CONCLUSIONS 

We have obtained experimental and numerical flux 
profiles that render the major features of flux flow at 
"vortex-lattice melting" in the presence of pinning and 
magnetic hysteresis. The flux profiles show a second front 
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corresponding to the well-known second peak in the mag- 
netization curve. In thin superconducting plates, "melt- 
ing" occurs well before the arrival of the second flux front. 
At that moment, the phase transition is complete and the 
low- field phase ( "vortex solid" ) with lower J c completely 
disappears there. The appearance of the high field phase 
( "vortex liquid" ) does not produce any sharp feature on 
the observable field profiles. Rather, there is a large re- 
gion preceding the second flux front in which domains 
of "solid" and "liquid" coexist dynamically: the relative 
width of these domains is controlled by the electrical field 
induced by the vortices moving under the force exerted 
by the current, which flows continuously across the do- 
mains; the current in its turn is controlled by the over- 
all electrodynamics of the system. This intricate feature 
should be taken into account in a correct interpretation 
of experiments. 

Evidently this effect occurs for any flux dynamics in 
the presence of a first order phase transition, irrespective 
of whether it is driven by ramping the applied field as in 
the presented experiment, or it happens during magnetic 
relaxation as in Ref. |f(|[ll| Moreover, when the process is 
slow and the current in the "weak" "solid vortex-lattice" 
phase is allowed to relax to zero, the "melting" proceeds 
nearly immediately everywhere because the flux flow re- 
quired for this is very low. In this case the remaining 
"solid" lamella should be very narrow and the second 
front moves on top of the nearly homogeneous "liquid" 
without any relation to the phase transition itself. 
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